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Abstract 

In this paper we study 0(2000) (quenched) lattice configurations from the APE 
collaboration, for different lattice volumes and for 6.0 < f3 < 6.4 using both the 
Wilson and the SW-Clover fermion actions. We determine the light hadronic spec- 
trum and meson decay constants and study the mesonic dispersion relation. We 
extract the hadronic variable J and t he s trange quark mass in the continuum at 
the next-to-leading order obtaining mf s (ii = 2GeV) = 122 ± 20 MeV. A study 
is made of their dependence on lattice spacing. We implement a newly developed 
technique to extract the inverse lattice spacing using data at the simulated values 
of the quark mass (i.e. at masses around the strange quark mass). 
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1 Introduction 



The lattice technique has proved a very effective theoretical tool to determine 
phenomenological quantities such as the mass spectrum and weak decay ma- 
trix elements. Unlike other approaches, it does not (in principle) suffer from 
uncontrolled approximations. However, in practice, one is forced to work on a 
lattice of (i) finite size, with (ii) a finite lattice spacing, and (iii) unphysically 
large masses for the light quarks. Also the quenched approximation is often 
used in lattice studies. In this work we determine the light hadronic spectrum 
and meson decay constants and we study the effects of (i), (ii) and (iii), and 
in particular, we try to overcome some of the inherent approximations cor- 
responding to (iii). It is important to study the systematics due to (i), (ii) 
and (iii) before the effects of the quenched approximation can be correctly 
understood. 

The results reported in this paper are obtained from (9(2000) (quenched) lat- 
tice configurations from the APE collaboration, for different lattice volumes 
and for 6.0 < (3 < 6.4 using both the Wilson action and the SW-Clover fermion 
action. In this range of j3 we cannot draw any conclusion about a dependence 
of hadron masses and pseudoscalar decay constants. On the other hand we 
find that the quenched approximation gives a resonable agreement with the 
experimental values. 

The plan of the paper is as follows. In the next section details of the simula- 
tions are given. Sect. 3 describes the method used to determine the spectrum 
and matrix elements from the lattice data. The results on hadron masses, me- 
son decay constants, J values, meson dispersion relations and quark masses 
are then presented in Sect. 4, with a comparison with results from other col- 
laborations. In the last section we report our concluding remarks. 



2 Lattice Details 

In the last four years the APE group has been extensively studying lattice 
QCD in the quenched approximation. The lattice data used in this study 
come from (9(2000) configurations generated by the APE collaboration pri- 
marily to study weak matrix elements such as /d, Jb and Bk and semileptonic 
decays [l]-[7]. The data have been obtained from eight sets of configurations 
with (3 = 6.0, 6.2 and 6.4. The quark propagators have been obtained from 
either the Wilson action or the 'improved' SW-Clover action, by the standard 
APE inversion procedure[8]. The hopping parameters and the other parame- 
ters used in each simulation are listed in Table 1. These K values correspond 
to quark masses roughly in the range of the strange quark mass which are 
adequate for the primary goal of the runs described before. In the past the 
APE collaboration has introduced the so called 'smearing' technique [8] in or- 
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Table 1 

Summary of the parameters of the runs analyzed in this work. T-Sweeps stands for 
thermalising Montecarlo sweeps performed starting from a cold configuration and 
I-Sweeps stands for the number of sweeps after which independent configurations 
are gathered. 





C60 


W60 


C62a 


W62a 


C62b 


W62b 


W64 


C64 


Ref 


[l]-[3] 


[1,4] 


[6] 


[6] 


[7] 


[4,5] 


[6] 


[6] 




6.0 


6.0 


6.2 


6.2 


6.2 


6.2 


6.4 


6.4 


Action 


sw 


Wil 


SW 


Wil 


SW 


Wil 


Wil 


SW 


# Confs 


200 


320 


250 


250 


200 


110 


400 


400 


Volume 


18 3 x 64 


18 3 x 64 


24 3 x 64 


24 3 x 64 


18 3 x 64 


24 3 x 64 


24 3 x 64 


24 3 x 64 


T-Sweeps 


3000 


3000 


5000 


5000 


3000 


3000 


7000 


7000 


I-Sweeps 


800 


800 


2000 


2000 


800 


2400 


3000 


3000 


K 






0.14144 


0.1510 






0.1488 


0.1400 




0.1425 


0.1530 


0.14184 


0.1515 


0.14144 


0.1510 


0.1492 


0.1403 




0.1432 


0.1540 


0.14224 


0.1520 


0.14190 


0.1520 


0.1496 


0.1406 




0.1440 


0.1550 


0.14264 


0.1526 


0.14244 


0.1526 


0.1500 


0.1409 



der to have a better determination of the lowest state mass. The APE group 
has compared data extracted from smeared and non-smeared propagators and 
found no real improvement on lattices of a time extension of 64 sites at f3 — 6.0 
and (3 = 6.2 [5]. The data presented in this work are obtained without any 
smearing procedure. In fact the data from W64 and C64 seem to suggest that 
a smearing procedure would have increased the reliability of the extraction of 
the mass of the lowest lying states. We plan to increase the time extension of 
this lattice and check the validity of the 'local' data presented here. 
A further comment is due regarding the 'thinning' procedure [8] used to reduce 
memory occupation. In this technique, correlation functions are computed tak- 
ing only one point out of three in each space direction. We have checked that 
this procedure does not affect the physical quantities extracted from the two 
point correlation functions. All baryonic data presented in this paper come 
from 'thinned' propagators while mesonic correlation functions come from the 
full lattice. 

Note that the simulations have been performed at (3 values of 6.0 or larger. 
This is to negate the large systematic error present in lattice data for (3 ^ 6.0 
due to lattice artifacts. A clear signal for these systematics is the lack of scal- 
ing amongst different physical quantities (see, for example, table 2 of [9] and 
[10]). 

Hadron masses and decay constants have been extracted from two-point cor- 
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relation functions in the standard way. For the meson masses and decay con- 
stants we have computed the following propagators: 



G55(t)=E( P 5(x,t)P 5 t (0,0)), 

x 

G 05 (t)=Y,(A (x,t)P^0,0)) , (1) 

X 

where 

P 5 (x,t) =iq(x,t)~f 5 q(x,t) , 
A li (x,t)=q(x,t)'y li 'y 5 q(x,t) . 

and the following propagators of the vector states: 

Gdt)= (2) 

i=l,3 x 

where 

Vi(x) = q(x,t)~fiq(x,t) . 

In order to determine the baryon masses we have evaluated the following 
propagators: 

G n (t)=52(N(x,t)Nl(0,0)) , 

X 

G s (t)=J2(A,(x,t) At (0,0)), (3) 

X 

where 

N = t abc {u a C^d b )u c 
A^ = e abc {u a C lfl u b )u c . 

For some meson correlation functions, non-zero momentum values were also 
studied (see section 4.2). 

We fit the zero- momentum correlation functions in eqs. 1, 2 and 3 to a sin- 
gle particle propagator with cosh or sink in the case of mesonic and axial- 
pseudoscalar correlation functions and with an exp function in the case of the 
baryonic correlation functions 



G 55 (t) = — exp(--M P sT) cosh(M P5 (- - *)) , 
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Table 2 

Time interval for baryons and mesons used in the fits for all lattices and for all 
momenta. 







C60 


W60 C62a W62a C62b W62b 


W64 


C64 








Mesons with zero momentum 










15-28 


15-28 18-28 18-28 18-28 18-28 


24-30 


24-30 








Mesons with non-zero momentum 






p 2 


= 1 


10-20 


10-20 12-20 12-20 12-20 


17-27 


17-27 


p 2 


= 2 


10-20 


10-20 11-16 11-16 10-15 


17-22 


17-22 


p 2 


= 3 


10-15 


10-15 11-16 11-16 












Baryons with zero momentum 










12-21 


12-21 18-28 18-28 18-28 18-28 


22-28 


22-28 



1 T 
Gu(t) = — exp(--M y T) cosh(M v (- - t)) , 

(4) 

G 05 (t) = — exp(--M PS T) sinh(M P5 (- - t)) , 

G nt s(t)=C n ' 5 exp(-M ntS t) , 

in the time intervals reported in Table 2. In eqs. 4, T represents the lattice time 
extension, the subscripts PS and V stand for pseudoscalar and vector meson, 
n and 5 stand for nucleon- and delta-like baryonsQ. To improve stability, 
the meson (axial-pseudoscalar) correlation functions have been symmetrized 
(anti-symmetrized) around t = T/2. The time fit intervals are chosen with 
the following criteria: we fix the lower limit of the interval as the one at 
which there is a stabilization of the effective mass, and, as the upper limit 
the furthest possible point before the error overwhelms the signal. We discard 
the possibility of fitting in a restricted region where a plateau is present, as 
the definition of such a region is questionable [11]. For lattices with highest 
number of configurations, i.e. W60, W64 and C64, we confirm that a higher 
statistics does not lead to a longer or better (relative to statistical errors) 
plateau [11]. The results of these fits for each run are reported in Tables 3-10. 
Representative examples of the effective mass plots can be seen in Figs. 1 - 4 in 
the case of W60 for the heaviest K = 0.1530. The errors have been estimated 
by a jacknife procedure, blocking the data in groups of 10 configurations and 
we have checked that there are no relevant changes in the error estimate by 
blocking groups of configurations of different size. We have also attempted to 

1 i.e. n stands for the nucleon N, AS or S baryon, and 5 is either a A ++ or a fi. 
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fit the first excited state mass for the mesons with a two hyperbolic cosine fit 
but the results are rather unstable and noisy on all lattices and will not be 
reported here. A similar conclusion holds for the baryonic parity partners. 

The pseudoscalar and vector decay constants f P $ and 1/ fy are defined through 
the equations 



(0\A \PS) = i^M PS , (5) 

A 

M 2 

(0\mr) = el—^, (6) 
Jv^v 

where e\ is the vector-meson polarization, M PS and My are the pseudoscalar 
and vector masses and Z v ,a are the renormalization constants. We extract fpg 
from the ratio 



Rfps(t) = ZA W) 

~^Z A — tanh(M P5 (--0) 

= ^Ptanh(M PS (|-()), (7) 

and the vector-meson decay constant is obtained directly from the parameters 
of the fit to Gu(t), eqs. 4: 



Zvfi 



v 



My 



(8) 



The results of these fits for each run are reported in Tables 3-10. 



3 Extraction of physical quantities 



Once the hadronic correlation functions have been fitted, and the lattice 
masses and matrix elements extracted, the 'conventional method' of determin- 
ing physical quantities relies on a number of chiral extrapolations. The first 
extrapolation is in the lattice pseudoscalar mass squared, (M PS a) 2 , against 
1/ K. This sets the critical kappa value K c corresponding to zero quark mass. 
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A second chiral extrapolation is performed of (Mya) to zero quark mass. Us- 
ing the experimental value of the p mass, My(K c ) = 770 MeV, the lattice 
spacing a is determined (assuming a massless pion). 

This conventional method relies on an extrapolation from a region where data 
has been obtained to a region well beyond the reach of direct simulations, 
assuming that the pseudoscalar mass follows the predictions of PCAC, and 
that other physical quantities scale linearly in the quark massQ- This problem 
suggests one should extract as much physics as possible from the 'strange' 
region and to use chiral extrapolations only where absolutely necessary. The 
method we use to achieve this aim is outlined below. 

• We define the following five 'lattice physical planes' for masses and decay 
constants: {My a, (M PS a) 2 ), (f PS a/Z A , (M PS a) 2 ), (l/(f v Z v ), (M PS a) 2 ), 
(M n a, (M PS a) 2 ) and (M s a, (M PS a) 2 ). We plot the lattice data for each K 
value used in the simulation on these planes . (See Figs. 5-9.) 

• On the (M v a , (M PS a) 2 ) plane we impose the physical ratio C s i = Mk* /Mk- 
This corresponds to find the intercept of the curve 



with the linear fit of the lattice data (i.e. the right-hand curve in Fig. 5). 
Lacock and Michael [13] have used this approach to define the quantity J. 
(See also [14] and [15].) The intercept of this curve and the simulated data 
defines M^a and Mx*a but only one of these is an independent prediction. 

• We now use the value of Mkcl determined above to read off the lattice meson 
decay constants, (/a-os/Z^) and (fx^Zy)' 1 from the corresponding fpg and 
fy planes. At this stage we have the following 3 independent predictions: 
(M K *a), (f K a/ZA) and (fx*Zy)~ 1 . These were all obtained from the quark 
mass region actually simulated (see Figs 5, 8 and 9), and hence they do not 
rely on a chiral extrapolation.^ 

• The lattice spacing a, can simply be set by comparing either (M^.a) or 
(Jkci/Za) with their experimental value, (/a* cannot be used since it is 
dimensionless.) Of these two quantities, we choose Mk* to set the scale 
whenever we need to quote results in physical units. This is because Mk* is 
determined experimentally more accurately than fx, and it has no ambigu- 
ities associated with a renormalisation constant. However, for completeness, 
we list in Table 11 the a" 1 values obtained from all the quantities. 

In this analysis for M K *a, /ku/Za an d (Jk'Zv)' 1 we have assumed that the 
quadratic 577(3) —flavour breaking effects in the meson masses or decay con- 

2 It is unlikely that the problem of 'chiral logs' in the quenched approximation be 
visible at the values of the simulated quark masses [12]. 

3 Depending on the choice of K, a modest linear extrapolation in M PS may be 
required. The range of this extrapolation is small and so should not cause problems. 




(9) 
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stants can be safely neglected in the range of masses explored here. This means 
that we assume that the masses and decay constants for mesons, with both 
degenerate, and non-degenerate constituent quarks, lie on the same universal 
curve. Studies have shown this to be an entirely reasonable hypothesis for 
quarks around the strange mass [16]. 

To proceed to obtain other physical quantities, we need the additional input 
of quark model considerations, and chiral extrapolations. Therefore the de- 
termination of these additional quantities suffers from the same criticisms as 
outlined above for the conventional approach. Consider the hadronic physical 
quantity M (where M = My, fp$, fy,M n or Ms). As a first approximation, 
the dependence of M on the hadron's quark masses can be approximated by 

Q 

Ma = A M + B M J2m q a, (10) 

q=l 

where the m q are the masses of the constituent quarks, and Q = 2, 3 for 
mesons, baryons. This is the usual assumption of the quark model. In order 
to continue we need the additional assumption of PCAC 

(M PS a) 2 = D(mia + m 2 a). (11) 

Combining eqs. 10 and 11 and eliminating the quark masses we have in the 
case of M = M v , 

M v a = A Mv + ^{M PS af. (12) 

Assuming this remains valid from the simulation region down to (near-) zero 
quark masses, we can obtain a prediction for M p a and (M w a) 2 as follows. 

• On the {Mya , (M PS a) 2 ) plane we impose the physical ratio Cu = M p /M w . 
The intercept of this curve and the linear extrapolation of the simulated 
data defines M p a and M w a (only one independent). 

• We now use the value of M n a determined above to read off the lattice meson 
decay constants, (/n-a/Z^) and (fpZy) -1 from the corresponding f P $ and 
fv planes. 

From Eq.(ll) we immediately have 

(M^a) 2 = 2(M K a) 2 - (M n a) 2 , (13) 

where i] ss is the hypothetical pseudoscalar meson with quark content ss as- 
suming no mixing with other states, for whose 'experimental' mass we have 
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M Vss = 0.686 GeV. Using the lattice value of M Vss a, we can determine M^ 
and (f^Zy)^ 1 from Figs. 5 & 9. As a check of the assumption of PCAC and 
Eq. 10, we plot in Fig. 10 the lattice values of M V /M K * against Mp S /M^ 
together with their experimental values. The data lie on a straight line to 
remarkable accuracy. 

For baryonic masses, the situation is more complicated. Eq.(10) can be written 
as: 

M n a = A Mn + B Mn (m 1 a + m 2 a + m 3 a), (14) 
M s a = A Ms + B M& {mia + m 2 a + m 3 a). 

Combining these with eq. (11) leads to 



M N a = 




+ 


D 


|(M.a) 2 






M As a = 




+ 


D 


{M K af + 






M s a = 




+ 


B M n 

D 


2(M K a) 2 






M A a = 


A Ms 


+ 


B Ms 
D 








M n a = 


A Ms 


+ 


B Ms 
D 


3(M K a) 2 







where 

M AS = ^M s + ^M A (15) 

The baryonic masses can be extracted by fitting the data for M n a and M$a 
against (M PS a) 2 to a straight line and using the values of M w a and M K a 
already determined, (see Figs. 6 and 7). Due to the limited number of K 
values, we do not attempt the fit to more complicated functions [17]. We have 
also checked that a quadratic fit does not give better x 2 - O ur predictions of 
the initial 3 quantities which do not require any chiral extrapolation have 
then been extended to include the following quantities which do require chiral 
extrapolations and quark model considerations: 

M p , M , U/Z A , {f p Z v )~\ {UZ v )-\ M N , M A2 , M s , M A andMf,. 

The method outlined above has the following properties: 

• We use only the lattice physical planes and not the unphysical planes where 
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the a;— axis is the quark mass 1/K. These unphysical planes will only be 
required when we need to determine unphysical quantities such as the quark 
masses. 

• This method allows us to fix the value of M K *a, and therefore cr 1 without 
the need of any chiral extrapolation. It also enables us to determine fx and 
fx* without any extrapolation. 

• With the additional assumption of PCAC (Eq.(ll)) and the quark model 
(Eq.(lO)) we can derive the rest of the meson spectrum and decay constants, 
and the baryon spectrum. 

• The conventional approach can be seen as the limit C« — > oo of this ap- 
proach. 

In the next section we use the above method to determine the physical quan- 
tities from the lattice simulations. 



4 Results 

This section briefly lists the results obtained using the method described in 
the previous section. Tables 3-10 contain the lattice data at each K value 
for M PS a, M v a, M n a, M$a, fpsa/Z A and (fyZy)^ 1 from the fits outlined is 
Section 2. Also listed are the K c values (obtained using Eq.(ll) with m q = 
1/2(1/ K q — l/K c ) and the chiral extrapolants using a linear fit with Mp s . 
All the mesonic and baryonic masses are also plotted in figs 10-12 and show 
a very good consistency, within the errors, among the various lattices. This 
shows that we observe a scaling behaviour in our data. 

4-1 Hadron masses and decay constants 

Using the method described in Sect. 3, the meson masses and decay constants 
(in lattice units) are displayed in Table 12 and 13. The inverse lattice spacing, 
a -1 , obtained from the experimental values is displayed in Table 11. 
Similarly the baryon masses, are displayed in Table 14. Again, for comparison, 
the a~ l values obtained using these masses are shown in Table 11. Note that 
all errors quoted in Tables 12, 13 and 14 are statistical only. In Tables 15 - 17, 
a listing is given of the lattice predictions of the mass spectrum and decay 
constants from C60 - C64. The scale was set from the K* mass. 

Comparing the clover meson masses at f3 = 6.2 from lattices C62a, C62b and 
the UKQCD results [16] in fig. 13, we find very good agreement with UKQCD 
for Lattice C62a. For C62b, while we find unexpectedly larger fluctuations in 
the comparison, we seem to get a somewhat lower mass; there may be some 
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residual finite size effects in the vector meson mass, even though one would 
rather expect a higher mass in this case. This problem may also be present 
for the (5 = 6.4 data for which the physical volume is the same as in C62b. 
Further investigation at larger lattice sizes is necessary to make the situation 
clearer. At (3 = 6.0 we can also compare our Wilson data with [25], [11] and 
[18], see Table 4 and fig. 14. The agreement is very good for the pseudoscalar 
data for which we have K values in common. Our data for the vector mass at 
the lightest quark mass is slightly higher. We do not have a clear explanation 
for this. We only note that if we had fitted our data in the t region in the 
interval 11 — 21 as done in fig. 2 of ref. [18] our value for the vector mass 
would have been My = 0.427(4) which is in much better agreement with their 
value. On the other hand, for the baryon data, we find very good agreement 
with the old APE [25] data, while we find slightly higher values and larger 
fluctuations when comparing with JLQCD [11] and LANL [18] data as shown 
in fig. 15. 

As can be seen from Tables 15 - 16, there is good consistency of the physical 
predictions between the different simulations to within statistical errors. Con- 
sidering the fact that the errors shown in Tables 15 - 16 are statistical only, 
and that these can often be underestimated [11], the agreement is remarkable. 

To compare the lattice decay constants with the experimental ones, we have 
used a 'boosted' one-loop form of the renormalisation constants [19] - [21]: 



Wilson action Za = 1 — 0.134g|^<7 
Z v = 1-0.174^ 
Clover action Z A = 1 — 0.0177g|^ 
Z v = 1-0.10^ 

where g*^ = 6/fe with [23] - [24] 

P m =<U P i aq >P + 0.15. (16) 



The results are shown in Table 17 and figs. 16 & 19. For the ratio fps/My in 
figs. 16-figs. 17, we see again a global consistency of our different data both 
for Wilson and the SW-Clover action. The UKQCD data [16], also reported 
in fig. 17, with the same renormalization decay constant Za disagree with ours 
for two of the three points where the comparison is possible. For the third 
point there is agreement but within much larger error bars. We do not see a 
strong a dependence for either actions comparing lattices at different (3 (C60, 
C62a and W60, W62a). This confirms statements made in Sect. 2 about the 
statistical presence of scaling in our data. The experimental points do lie quite 
well in the band extrapolated/interpolated from the lattice data. 
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The values of the pseudoscalar decay constants in Table 17 show larger de- 
viations from the experimental data than the ratio fps/My. It is clear that 
a cancellation of systematic error occurs in this ratio. Still there is an agree- 
ment within 1.5 standard deviations with the experimental data apart from 
the W60 value which is quite high. In particular it is higher than the value 
obtained for the Wilson action at (3 = 6 in ref. [22]. There the authors use 
different renormalization schemes than the one used here but the discrepancy 
is too large to be attributed to this. 

The situation is more delicate for the vector decay constant for which a depen- 
dence on the volume and a~ l may be present but with our data it would be 
difficult to disentangle the two effects. In particular in figs. 18, 19 and Table 17 
the values for C62a and C62b (SW-Clover on different physical volumes), C62a 
and W62a (SW-Clover and Wilson with the same physical volume), W60 and 
W62a ( Wilson with the same physical volume but different values of (3) and 
C60 and C62a (SW-Clover on the same volume and different j3) seem to show 
this problem. The extrapolated values for 1/ fv suffer very much due to these 
effects. The 'strange' vector decay constants seem to be slightly more stable 
also because no extrapolation is needed. 



4-2 Meson Dispersion relations 



For all the lattices except for W62b the non-zero momentum meson correlation 
functions were also computed. We can study discretization errors and finite 
volume effects testing which dispersion relation is best verified from data. We 
have tested the following dispersion relations 



E 2 =p 2 + M 2 , 

sinh 2 ^) = sin 2 (^) + sinh 2 ( — ) , 
v 4 4 4 

sinh 2 (|)=sin 2 (|)+sinh 2 (f), 

sinh 2 (E) = sin 2 (p) + sinh 2 M , 



(17) 
(18) 

(19) 
(20) 



where (17) is the continuum relation, while the others are lattice dispersion 
relations from different choices of lattice action. M — E(p — 0) is the mass. 
From the lattice we have extracted the energy E(p) of the states fitting the 
non-zero momentum correlation functions with a cosh function. The results we 
obtained for E(p) for pseudoscalar and vector mesons are reported in Table 18. 
We see only small differences among the various dispersion relations for the 
quark masses we have simulated, but eq. (19) seems to be favoured [2]. This is 
consistent with the conclusion of ref. [18] which refers however to much larger 
quark masses. 
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4-3 Strange Quark Mass 



Lattice QCD is in principle able to predict the mass of any quark by fixing 
to its experimental value the mass of a hadron containing a quark with the 
same flavour. The 'bare' lattice quark mass m(a) can be extracted directly 
from lattice simulations and can be related to the continuum mass m MS (fj,) 
renormalized in the minimal-subtraction dimensional scheme through a well- 
defined perturbative procedure [26] - [28] . Following ref. [29] 

m m Qi) = Z^a)m(a) 

where m(a) is the bare lattice quark mass related to the hopping parameter 

as 

m(a)a = l(l - h ■ (21) 



and Z^ s (fia) is the mass renormalization constant at scale fi which we choose 
to be 2 GeV. The hopping parameters and the strange quark masses reported 
in Tables 19 and 20 are extracted using the eq. (21) and (11) from the lattice 
values M w a and M K a. The error on m MS (fi) is estimated as in ref. [29] tak- 
ing into account the spread due to different definitions of the strong coupling 
constant. This leads to somewhat larger errors for the Wilson data than for 
the Clover ones due to the smaller fluctuation of the different renormalization 
constants in the improved case. There is a good consistency among the val- 
ues coming from the different lattices. Fig. 20 shows that we do not see any 
dependence on a within the errors in the Wilson data and a mild tendency in 
the clover data to decrease with increasing f3. We therefore conclude that any 
0(a) effects present are beneath the level of statistics, and/or hidden among 
finite volume effects. We then extract the average value of the strange quark 
mass without any extrapolation and get: 

mf^in = 2GeV) = 122 ± 20Me^ 

which is in agreement with the result of ref. [29]. It is also compatible with 
the value mf^(/i = 2GeV) = 90 ± 20 MeV of ref. [30], but one should 
take into account that this value comes from an analysis on various Wilson 
and Staggered lattices at different values of (3 and an extrapolation in a. The 
quark masses can also be extracted from the Ward identity for the axial vector 
current. A full analysis of the light quarks masses from both methods will be 
presented in a forthcoming paper. 
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4-4 J value(s) 



In this section we analyze the effect of the quenched approximation using the 

J variables first proposed for mesons by Lacock and Michael [13]. 

In the spirit of the Maiani-Martinelli proposal [14] [16] to use (Mk*—M p ) / [M\— 

M%) as an independent way to determine the lattice scale and avoiding the 

chiral extrapolation as described in section 3, Lacock and Michael proposed 

that 



dM v 

L PS 



J = M 'W- < 22) 



be computed on the lattice and compared with experimental values. J is a 
combination of parameters which is independent of a and K and which does 
not require any chiral extrapolation to the chiral limit. We use the exper- 
imental values of the meson masses reported in Table 15 to determine the 
experimental values of J. This gives the dimensionless value reported in Ta- 
ble 21. 

Using the method described in section 3 it is very easy to extract the J value 
for mesons. The results we have obtained are reported in the Table 21. The 
lattice values for J are significantly different to the experimental value and 
one can see a dependence on j3 but as discussed before a volume effect may 
be present on C62b, W64 and C64 which prevents a clear conclusion. The dis- 
agreement between the data and the experimental value of J is also evident in 
fig. 10. It is easy to see that, in the plane of fig. 10, the J value is proportional 
to the slope of the line through the experimental data which is clearly different 
to the slope of the lattice data. 

As well as defining the quantity J for the meson sector, a corresponding quan- 
tity for the nucleon and delta sector can be introduced. This can be viewed 
as the slope of the data in figs. 11, 12. While the data is noisy, it is likely that 
the lattice data also underestimates the slope of the experimental points in 
these plots. 



5 Conclusions 



We have analyzed a large set of data on lattices of different lattice spacing and 
lattice sizes with both the Wilson and SW-Clover actions. These simulations 
were not originally performed for a hadron spectrum study, hence the rela- 
tively large values for the light quark masses. To overcome potential problems 
with the chiral extrapolation, a new technique of obtaining lattice predictions 
without any extrapolation in the quark mass was developed. 
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In the (3 range studied, we conclude that there is no statistical evidence for an 
a dependence in the hadron masses or for the pseudoscalar decay constants. 
Our data agrees with the experimental data to within ~ 5% for mesonic 
masses and ~ 10% — 15% for baryonic masses and the pseudoscalar decay 
constants. A larger deviation is present for the vector decay constants which, 
in our opinion, deserve a more careful study at larger volumes and (3. The 
effects of finite size, if present, are within the statistical errors of our data. 

We obtain for the strange quark mass mf^(/i = 2GeV) = 122 ± 20 MeV 
without seeing any clear a dependence. The meson dispersion relations and 
the quantity J were also studied. 

We believe that the lattice results at (3 > 6 are quite stable and show a scaling 
behaviour, within current statistical errors. 
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Table 3 

Hadron masses and meson decay constants versus the hopping parameter for Lattice 



C60. 



K 


M PS a 


My a 


M n a 


M s a f PS a/Z A 


VtfvZv) 


0.1425 


0.439(1) 


0.549(5) 


0.833(9) 


0.90(2) 0.095(2) 


0.324(9) 


0.1432 


0.382(1) 


0.512(7) 


0.76(1) 


0.84(3) 0.089(2) 


0.34(1) 


0.1440 


0.308(1) 


0.48(1) 


0.69(1) 


0.78(5) 0.082(3) 


0.38(2) 


0.14549(2) 




0.40(2) 


0.54(2) 


0.65(7) 0.069(4) 


0.43(4) 


Table 4 

As in Table 3 for Lattice W60. 


K 


M PS a 


My a 


M n a 


Msa fpsa/ZA 


VifvZy) 


0.1530 


0.423(1) 


0.508(3) 


0.801(6) 


0.864(8) 0.114(2) 


0.391(6) 


0.1540 


0.364(1) 


0.468(4) 


0.729(7) 


0.81(1) 0.108(2) 


0.420(8) 


0.1550 


0.298(1) 


0.431(6) 


0.66(1) 


0.75(2) 0.100(2) 


0.46(1) 


0.15703(2) 




0.353(9) 


0.52(1) 


0.64(3) 0.087(3) 


0.52(2) 


Table 5 

As in Table 3 for Lattice C62a. 


K 


M PS a 


My a 


M n a 


M s a f PS a/Z A 


if(fvZv) 


0.14144 


0.297(1) 


0.386(3) 


0.586(5) 


0.066(1) 


0.310(6) 


0.14184 


0.258(1) 


0.362(5) 


0.542(8) 


0.062(2) 


0.326(9) 


0.14224 


0.215(1) 


0.34(1) 


0.50(1) 


0.057(2) 


0.35(2) 


0.14264 


0.162(2) 


0.34(3) 


0.46(4) 


0.052(3) 


0.40(6) 


0.14315(1) 




0.29(2) 


0.40(3) 


0.046(3) 


0.39(3) 


Table 6 

As in Table 3 for Lattice W62a. 


K 


M PS a 


My a 


M n a 


M s a f PS a/Z A 


1/ifvZy) 


0.1510 


0.289(1) 


0.366(2) 


0.566(4) 


0.080(1) 


0.381(6) 


0.1515 


0.254(1) 


0.343(3) 


0.525(6) 


0.075(1) 


0.401(7) 


0.1520 


0.215(1) 


0.321(5) 


0.48(1) 


0.070(2) 


0.42(1) 


0.1526 


0.158(1) 


0.29(1) 


0.45(3) 


0.062(2) 


0.45(3) 


0.15329(1) 




0.26(1) 


0.38(2) 


0.056(2) 


0.48(2) 
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Table 7 

As in Table 3 for Lattice C62b. 



K 


M PS a 


My a 


M n a 


M s a 


f PS a/Z A 


V(fvZy) 


0.14144 


0.294(3) 


0.378(6) 


0.59(2) 


0.68(1) 


0.066(2) 


0.29(1) 


0.14190 


0.249(3) 


0.34(1) 


0.53(2) 


0.65(2) 


0.061(2) 


0.30(2) 


0.14244 


0.186(4) 


0.30(3) 


0.43(5) 


0.57(4) 


0.055(3) 


0.29(4) 


0.14312(4) 




0.25(3) 


0.35(7) 


0.53(6) 


0.048(4) 


0.30(4) 


Table 8 

As in Table 3 for Lattice W62b. 


K 


M PS a 


My a 


M n a 


M s a 


fpsa/ZA 


V(fvZv) 


0.1510 


0.291(1) 


0.370(3) 


0.576(7) 


0.618(7) 


0.079(2) 


0.390(9) 


0.1520 


0.217(1) 


0.323(7) 


0.50(2) 


0.56(2) 


0.069(2) 


0.43(2) 


0.1526 


0.160(2) 


0.29(1) 


0.42(6) 


0.55(4) 


0.062(2) 


0.46(3) 


0.15330(2) 




0.26(1) 


0.38(6) 


0.49(4) 


0.055(2) 


0.49(3) 


Table 9 

As in Table 3 for Lattice W64. 


K 


M PS a 


My a 


M n a 


M s a 


f PS a/Z A 


1/ifvZv) 


0.1488 


0.235(1) 


0.289(2) 


0.461(4) 


0.499(7) 


0.0609(9) 


0.333(4) 


0.1492 


0.205(2) 


0.266(3) 


0.423(5) 


0.47(1) 


0.0573(9) 


0.346(6) 


0.1496 


0.172(2) 


0.242(4) 


0.38(1) 


0.42(2) 


0.053(1) 


0.359(8) 


0.1500 


0.136(3) 


0.22(1) 


0.33(3) 


0.36(4) 


0.048(2) 


0.38(2) 


0.15058(3) 




0.188(9) 


0.28(2) 


0.33(4) 


0.043(2) 


0.39(2) 


Table 10 

As in Table 3 for Lattice C64. 


K 


M PS a 


My a 


M n a 


M 5 a 


f PS a/Z A 


l/{fvZv) 


0.1400 


0.249(1) 


0.306(2) 


0.484(4) 


0.52(1) 


0.0524(9) 


0.269(5) 


0.1403 


0.220(1) 


0.283(3) 


0.446(6) 


0.49(1) 


0.049(1) 


0.276(6) 


0.1406 


0.188(2) 


0.259(5) 


0.40(1) 


0.44(2) 


0.046(1) 


0.281(8) 


0.1409 


0.152(3) 


0.23(1) 


0.34(2) 


0.37(4) 


0.042(1) 


0.29(2) 


0.14143(3) 




0.19(1) 


0.29(2) 


0.32(5) 


0.037(2) 


0.30(2) 
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Table 11 

Values of a -1 in GeV from different observables for all lattices. The values of a~^\ 
(shown in boldface) are taken as our preferred estimate of the inverse lattice spacing. 





C60 


W60 


C62a 


W62a 


C62b 


W62b 


W64 


C64 


a' 1 


1.88(9) 


2.15(5) 


2.6(2) 


2 9(1) 


3.0(3) 


2.9(2) 


4 0(2) 


3.9(2) 


(1, 


1.98(8) 


2.26(5) 


2.7(1) 


3.00(9) 


3.0(3) 


3.0(1) 


4.1(2) 


4.0(2) 




1.98(8) 


2 26(5) 


2 7(1) 


3 00(9) 


3.0(3) 


3 0(1) 


4 2(2) 


4 0(2) 


a 4> 


2 06(7) 


2 35(5) 


2 8(1 ) 


3 1 0(8) 

tj. -L\J I O I 


3 1 (3) 


3 1(1) 


4 3(2) 


4 1 (2) 


a N~ 


1.70(6) 


1.79(5) 


2.3(1) 


2.4(1) 


2.7(5) 


2.4(3) 


3.2(2) 


3.2(2) 


a AS 


1.86(6) 


1.99(4) 


2.5(1) 


2.7(1) 


2.8(4) 


2.7(3) 


3.6(2) 


3.5(2) 




1.98(6) 


2.13(4) 


2.68(8) 


2.84(8) 


2.9(4) 


2.8(2) 


3.8(2) 


3.7(2) 




1.9(2) 


1.90(8) 






2.3(2) 


2.5(2) 


3.7(4) 


3.8(5) 




2.07(9) 


2.21(6) 






2.7(1) 


2.9(1) 


4.0(2) 


4.0(3) 




1.9(1) 


1.90(7) 


2.9(2) 


2.9(1) 


2.8(2) 


2.9(1) 


3.7(2) 


3.6(2) 


IK 


2.11(8) 


2.16(6) 


3.0(1) 


3.1(1) 


3.1(2) 


3.17(9) 


4.1(1) 


4.0(1) 



Table 12 

Extrapolated/interpolated meson masses in lattice units for all lattices. 





M p a 


Mx*a 


M Vss a 


M a 


C60 


0.41(2) 


0.45(2) 


0.35(1) 


0.49(2) 


W60 


0.358(9) 


0.395(9) 


0.303(7) 


0.433(9) 


C62a 


0.30(2) 


0.33(2) 


0.25(1) 


0.36(1) 


W62a 


0.27(1) 


0.298(9) 


0.229(7) 


0.329(9) 


C62b 


0.26(3) 


0.29(3) 


0.23(2) 


0.33(3) 


W62b 


0.26(1) 


0.30(1) 


0.23(1) 


0.33(1) 


W64 


0.192(8) 


0.215(8) 


0.165(6) 


0.239(8) 


C64 


0.20(1) 


0.22(1) 


0.171(8) 


0.25(1) 
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Table 13 

Extrapolated/interpolated meson decay constants for all lattices. 

f n a/(Z A m p a) l/(f p Z v ) f K a/{Z A m K *a) \/(f K *Z v ) l/(f (j) Z v ) fx/U 



C60 


0.17(1) 


0.42(3) 


0.172(9) 


0.39(2) 


0.36(1) 


1.11(2 


W60 


0.25(1) 


0.51(2) 


0.239(8) 


0.48(2) 


0.45(1) 


1.08(1 


C62a 


0.16(1) 


0.39(3) 


0.164(9) 


0.36(2) 


0.332(8) 


1.15(3 


W62a 


0.21(1) 


0.47(2) 


0.214(8) 


0.45(1) 


0.417(8) 


1.13(1 


C62b 


0.19(3) 


0.30(4) 


0.18(2) 


0.30(3) 


0.30(2) 


1.10(3 


W62b 


0.21(1) 


0.49(3) 


0.21(1) 


0.46(2) 


0.43(1) 


1.12(2 


W64 


0.23(2) 


0.39(2) 


0.22(1) 


0.38(1) 


0.363(8) 


1.10(2 


C64 


0.19(1) 


0.30(2) 


0.18(1) 


0.29(1) 


0.284(9) 


1.10(1 



Table 14 



Extrapolated/interpolated baryon masses in lattice units for all lattices. 





M N a 


M AS a 


M s a 


M A a 


M n a 


C60 


0.55(2) 


0.61(2) 


0.67(2) 


0.66(6) 


0.81(3) 


W60 


0.53(1) 


0.57(1) 


0.62(1) 


0.65(3) 


0.76(2) 


C62a 


0.41(3) 


0.45(2) 


0.49(1) 






W62a 


0.39(2) 


0.43(2) 


0.46(1) 






C62b 


0.35(7) 


0.40(6) 


0.45(5) 


0.53(6) 


0.62(2) 


W62b 


0.39(5) 


0.43(5) 


0.46(4) 


0.50(4) 


0.57(2) 


W64 


0.29(2) 


0.32(2) 


0.35(1) 


0.34(4) 


0.42(2) 


C64 


0.29(2) 


0.32(2) 


0.35(2) 


0.33(5) 


0.42(3) 
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Table 15 

Predicted meson masses in GeV for all lattices, using the scale from Mk* ■ 





M p 


M Vaa 




Exper. 


0.770 


0.686 (see text) 


1.019 


C60 


0.809(7) 


0.6849(3) 


0.977(7) 


W60 


0.808(3) 


0.6849(1) 


0.978(3) 


C62a 


0.81(1) 


0.6849(5) 


0.98(1) 


W62a 


0.803(6) 


0.6851(2) 


0.984(6) 


C62b 


0.79(1) 


0.6856(5) 


1.00(1) 


W62b 


0.797(7) 


0.6853(3) 


0.989(7) 


W64 


0.796(4) 


0.6853(2) 


0.990(4) 


C64 


0.792(4) 


0.6855(2) 


0.994(4) 



Table 16 

Predicted baryon masses in GeV for all lattices, using the scale from Mk* ■ 





M N 


Mas 


M s 


M A 


M n 


Exper. 


0.9389 


1.135 


1.3181 


1.232 


1.6724 


C60 


1.09(5) 


1.21(4) 


1.32(4) 


1.3(1) 


1.60(9) 


W60 


1.19(5) 


1.29(4) 


1.40(4) 


1.46(7) 


1.71(4) 


C62a 


1.1(1) 


1.22(8) 


1.34(7) 






W62a 


1.17(7) 


1.28(6) 


1.39(5) 






C62b 


1.1(2) 


1.2(2) 


1.4(1) 


1.6(3) 


1.9(2) 


W62b 


1.2(1) 


1.3(1) 


1.40(9) 


1.50(9) 


1.72(5) 


W64 


1.21(9) 


1.32(8) 


1.43(6) 


1.4(2) 


1.72(9) 


C64 


1.2(1) 


1.29(8) 


1.41(7) 


1.3(2) 


1.7(1) 
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Table 17 



Predicted meson decay constants for all lattices, using the scale from Mr-*. 





U (GeV) 


1 

Tp 


Sk {GeV) 


1 

fx* 


1 

u 


Exper. 


0.1307 


0.28 


0.1598 




0.23 


C60 


0.134(9) 


0.35(3) 


0.149(8) 


0.33(2) 


0.30(1) 


W60 


0.155(7) 


0.37(2) 


0.167(6) 


0.35(1) 


0.324(8) 


C62a 


0.124(9) 


0.33(3) 


0.143(8) 


0.30(2) 


0.281(7) 


W62a 


0.135(6) 


0.35(2) 


0.153(5) 


0.33(1) 


0.307(6) 


C62b 


0.14(2) 


0.26(4) 


0.16(2) 


0.25(3) 


0.25(2) 


W62b 


0.135(8) 


0.36(2) 


0.152(7) 


0.34(2) 


0.315(9) 


W64 


0.147(9) 


0.29(1) 


0.161(8) 


0.283(9) 


0.272(6) 


C64 


0.144(9) 


0.25(1) 


0.158(8) 


0.25(1) 


0.243(8) 
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Table 18 



Values of E(p) for mesons with non-zero momentum versus the hopping parameter. 



K p 2 = 1 p 2 = 2 p 2 = 3 


p 2 = 1 p 2 = 2 p 2 = 3 


C60 0.1425 0.567(4) 0.67(1) 0.73(3) 
0.1432 0.526(6) 0.63(2) 0.68(4) 
0.1440 0.479(9) 0.57(4) 0.60(8) 


0.654(4) 0.751(8) 0.82(1) 
0.625(6) 0.73(1) 0.78(1) 
0.60(1) 0.70(2) 0.70(2) 


W60 0.1530 0.554(5) 0.66(1) 0.73(2) 
0.1540 0.513(7) 0.62(2) 0.68(4) 
0.1550 0.47(1) 0.58(3) 0.62(7) 


0.621(4) 0.716(7) 0.79(1) 
0.591(6) 0.692(9) 0.75(2) 
0.566(8) 0.67(2) 0.69(3) 


C62a 0.14144 0.399(3) 0.479(6) 0.53(2) 
0.14184 0.371(4) 0.455(9) 0.51(4) 
0.14224 0.340(7) 0.43(1) 0.50(7) 
0.14264 0.30(1) 0.40(3) 


0.476(3) 0.557(4) 0.618(8) 
0.455(3) 0.543(7) 0.60(1) 
0.436(5) 0.54(1) 0.59(2) 
0.42(1) 0.56(3) 0.55(3) 


W62a 0.1510 0.395(3) 0.479(6) 0.52(2) 
0.1515 0.369(3) 0.459(9) 0.50(3) 
0.1520 0.343(5) 0.44(1) 0.48(5) 
0.1526 0.31(1) 0.41(3) 0.5(1) 


0.459(2) 0.541(3) 0.605(7) 
0.440(3) 0.528(5) 0.592(9) 
0.423(3) 0.518(8) 0.58(1) 
0.404(6) 0.53(2) 0.54(3) 


C62b 0.14144 0.48(1) 0.60(3) 
0.14190 0.46(2) 0.56(4) 
0.14244 0.42(6) 0.48(7) 


0.536(7) 0.66(2) 0.77(4) 
0.53(1) 0.64(3) 0.77(5) 
0.53(3) 0.63(5) 0.76(8) 


W62b - 




W64 0.1488 0.363(6) 0.47(2) 
0.1492 0.35(1) 0.46(4) 
0.1496 0.33(2) 0.46(8) 
0.1500 0.33(5) 0.4(2) 


0.406(4) 0.50(1) 
0.394(5) 0.50(2) 
0.385(8) 0.51(4) 
0.38(2) 0.6(1) 


C64 0.1400 0.371(6) 0.49(2) 
0.1403 0.36(1) 0.49(4) 
0.1406 0.34(2) 0.50(8) 
0.1409 0.33(4) 0.5(2) 


0.419(4) 0.51(1) 
0.406(5) 0.51(2) 
0.393(8) 0.52(5) 
0.38(1) 0.6(2) 
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Table 19 



Values for K c , 


K\ and K s 


for all lattices. 






K c 


Ki 


K s 


C60 


0.14549(2) 


0.14540(2) 


0.1436(1) 


W60 


0.15703(2) 


0.15693(2) 


0.15492(9) 


C62a 


0.14315(1) 


0.14310(2) 


0.1419(1) 


W62a 


0.15329(1) 


0.15322(2) 


0.15184(9) 


C62b 


0.14312(4) 


0.14308(4) 


0.1421(2) 


W62b 


0.15330(2) 


0.15323(2) 


0.1519(1) 


W64 


0.15058(3) 


0.15054(3) 


0.14969(8) 


C64 


0.14143(3) 


0.14140(3) 


0.14074(7) 



Table 20 

Values for the lattice strange lattice quark masses for all lattices and the corre- 
sponding MS values at NLO, both in MeV. 





m s (a) 


mf s {n = 2 GeV) 


C60 


89(3) 


120(10) 


W60 


98(2) 


130(20) 


C62a 


83(4) 


120(10) 


W62a 


93(3) 


130(10) 


C62b 


75(6) 


110(10) 


W62b 


92(4) 


130(20) 


W64 


82(3) 


120(10) 


C64 


69(3) 


100(10) 
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Table 21 

Predicted values of J for all lattices. 



J 


Exper. 


0.48 


C60 


0.34(3) 


W60 


0.34(1) 


C62a 


0.34(5) 


W62a 


0.36(2) 


C62b 


0.42(5) 


W62b 


0.38(3) 


W64 


0.39(2) 


C64 


0.40(2) 
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Fig. 1. Effective pseudoscalar mass Mpg in lattice units versus time for lattice W60 

and K = 0.1530. Dashed lines show the fitted value and its error in the chosen time 

interval. 
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Fig. 2. As in fig. 1 for the vector meson mass My. 



26 



0.900 



0.875 



0.850 



0.825 — 



0.800 



0.775 



n 1 1 1 1 1 1 r 







* 






T > 









X >C 



750 ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 

10 15 20 25 30 

t/a 

Fig. 3. As in fig. 1 for the nucleon mass M n . 
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Fig. 4. As in fig. 1 for delta mass Mg. 
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Fig. 5. Values of My versus Mp s in lattice units and the fitting line for W60. The 
curves represent the equations My a = Cy / (Mpsa) 2 with C = Cu (leftmost) and 
C = C s i (rightmost) defined in section 3. 
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Fig. 6. Values of M n versus Mp S in lattice units and the fitting line for W60. 
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Fig. 9. As in fig. 6 for 
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Fig. 10. Mass of vector meson My versus the squared pseudoscalar mass Mp S for ali 
iattices. My and Mps normalized to the lattice value of Mk* and Mk respectively. 
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Fig. 11. As in fig. 10 for M n . 
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Fig. 12. As in fig. 10 for M s . 
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Fig. 13. Vector meson masses in lattice units at f3 = 6.2 for SW-Clover action from 
UKQCD [16], Lattice C62a and Lattice C62b. 
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Fig. 14. Vector meson masses in lattice units at (3 = 6.0 for Wilson action from 
Lattice W60, APE 24 3 x 32 [25], JLQCD [11] and LANL [18]. 
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Fig. 15. Nucleoli masses in lattice units at (3 = 6.0 for wilson action from Lattice 
W60, APE 24 3 x 32 [25], JLQCD [11] and LANL [18]. 
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Fig. 16. Values of the ratio versus the squared pseudoscalar mass Mp S nor- 
malized to the lattice value of Mk for all lattices with Wilson action. 
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Fig. 17. Values of the ratio versus the squared pseudoscalar mass Mp s nor- 
malized to the lattice value of Mk for all lattices with SW-Clover action. 
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Fig. 18. Values of the vector decay constant 1/fv versus the squared pseudoscalar 
mass Mp S normalized to the lattice value of Mk for all lattices with Wilson action. 
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Fig. 19. Values of the vector decay constant 1/ fv versus the squared pseudoscalar 
mass Mp S normalized to the lattice value of Mk for all lattices with SW-Clover 
action. 



0.16 



0.14 — 



0.12 — 



0.10 — 



0.08 



060 o 
W60 □ 
062a + 
W62a w 
062b * 
W62bx 
W64 ;:s 
064 o 



0.1 0.2 0.3 

a [GeV -1 ] 



0.4 



0.5 



0.6 



Fig. 20. Renormalized quark masses m^f s for all lattices. 
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